Quantum oscillation in the Hofstadter Model 
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We have developed a novel quantum transfer matrix method to study thermodynamic properties 
of the Hofstadter model on square lattices. This method overcomes the finite size effect and allows 
the magnetization and other thermodynamic quantities to be accurately and efficiently calculated. 
The interplay between the lattice potential and the applied magnetic field leads to substantial 
corrections to the conventional de Haas-van Alphen (dHvA) oscillation of magnetization. It yields 
many paramagnetic peaks or diamagnetic deeps above the background of dHvA oscillation. This 
correction bears the fractal feature of energy spectra and can be used to probe the hierarchical 
structure of Hofstadter butterfly by experimental measurements. 



PACS numbers: 71.18.-fy, 02.30.1k, 75.20.-g 

The interplay between a uniform magnetic field and 
a periodic crystal potential on two-dimensional electron 
gas has attracted the condensed matter community for 
more than three decades [l], H, H, 0, H)]. Great atten- 
tion has been paid to the investigation of the Hofstadter 
model[B], IH, in which the periodic potential is mimicked 
by tight binding electrons. This model has been solved by 
the Bethe Ansatz0,|I| and exact diagonalizationQ meth- 
ods when the magnetic flux quanta <f> penetrating through 
a unit cell is rational <f> = p/q (p and q are mutually prime 
integers) with relatively small q. The band structure is 
found to have a fractal butterfly shape. This topological 
structure of spectra, known as the Hofstadter butterfly, 
was observed in the Hall conductance measurement of 
two-dimensional electron gas in a superlattice [Toj . How- 
ever, to study thermodynamic properties, especially the 
lattice correction to the quantum oscillation of electrons 
in a magnetic field, one has to solve this model for ar- 
bitrary <p. This is a very challenging problem since the 
largest q that can be handled by the Bethe Ansatz or 
exact diagonalization is generally less than 1000. 

In the absence of periodic potential, the orbital magne- 
tization of an electron gas will oscillate with an applied 
magnetic field. This well known dHvA oscillation has 
been extensively used to detect Fermi surface volumes 
of metallic crystals. At zero temperature, the ground 
energy of the Hofstadter model is not an analytic func- 
tion of applied magnetic field due to the fractal feature 
of spectra. Thus the orbital magnetic susceptibility is 
not a well defined quantity at zero temperature. How- 
ever, thermal fluctuation can remove this non-analytical 
feature, and the orbital magnetization and susceptibility 
are still physically well defined quantities at finite tem- 
perature. 

In this Letter, we propose a novel quantum transfer 
matrix method to study the lattice correction to the 
dHvA oscillation in the Hofstadter model on square lat- 
tices at finite temperature. In the framework of Landau 
gauge, we will show that the partition function of the 
system can be expressed as a product of a series of 2 x 2 



local transfer matrices. This method avoids direct diag- 
onalization of the Hamiltonian and allows the thermo- 
dynamic limit to be explored directly and reliably. At 
low field and high temperature, the lattice correction is 
small, we find that the magnetization oscillates with ap- 
plied magnetic field, similar as in the system without 
periodic potentials. However, at high field and low tem- 
perature, the lattice correction becomes significant. In 
addition to the dHvA oscillation, many sharp peaks or 
deeps are present in both the magnetization and suscep- 
tibility curves. These peaks or deeps exhibit hierarchical 
structures. They result from the fractal band dispersion 
and can be used to probe the topological structure of 
Hofstadter butterfly in experiments. 

The Hofstadter model describes dynamics of tight 
binding electrons in a uniform magnetic field. In the Lan- 
dau gauge, if we assume the vector potential is nonzero 
only on bonds along the y-direction, the system will be 
translationally invariant along the x-direction. By tak- 
ing the plane wave expansion (or Fourier transformation) 
along the x-axis, one can then decouple the Hofstadter 
model H into a sum of a series of one-dimensional Hamil- 
tonian Hk 



H = $> fe 
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where k — 2nn/N x (n = 0, 1, ...N x — 1) is the momen- 
tum of electrons along the x-axis and N x is the lattice 
dimension along that direction, y is the lattice coordi- 
nate of electrons along the y-axis. <f> is the magnetic flux 
penetrating each plaquette. 

Given k, the partition function of Hk is defined by 

Z k = Trcxp(-/3i^), (3) 

where = 1/kgT and T is temperature. The partition 
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function of the whole system is simply a product of Zk 
for all k. 

To evaluate Z k , let us first divide H k into two parts, 

Hk = Hk.ev 



H k ,odd, where 

Hk,odd 



y 

^2 hk - 2 y 



and 

hk,y = t{ c l,y+i c k,y + h.c) + 2t cos (27ry0 - ft) 4, v c fc.i/- 

The individual terms in H k ^ even or Hk t odd commute 
with each other. Thus it is relatively simple to evaluate 
thermodynamic quantities of Hk, even or Hk.odd- To uti- 
lize this property, let us divide f3 into M equivalent parts, 
e — /3/M and apply the Trotter-Suzuki formula [III. Il2| 
to decompose Z k as 

Z k = Tr ( e - eH *.odd e -eH k „ 



0(e 2 



(4) 



afe.a = -et[cos (27r?/^ - ft) - /x/2], ^ k , y z 
u k , y = -etswhjk, y exp(ak,y)/lk,y, 



and 



\Jpk,y<lk,y + r k y - P" i s the chemical potential. In general, 
for any quadratic Hamiltonian, it can be shown that the 
corresponding local transfer matrix can be always writ- 
ten as an exponent of a quadratic function of fermion 
operators [13|. 

Reversing the order of I and y in Eq. ([5]), one can then 
reexpress the partition function as a product of transfer 
matrices 



Z k = lim Tr(ri i2 T 2 .3 ■ • ■ T N1 ), 
where T VtV+ i are transfer operators defined by 



(9) 
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By inserting completeness identities to the above expres- 
sion, we have 



M N y / 2 



where 



Z k 



1,1+1 



an n 

i=i «=i 



21-1,21 21,21+1 
V 2y-l,2y V 2y,2y+n 



(nln'\e~ eh ^ 



nf+V +1 



(5) 



(6) 



J y,y+1 ~ \ IL y IL y+i\^ l'"y '"y+ 

The subscripts represent the lattice positions and the su- 
perscripts represent the coordinates in the inverse tem- 
perature or Trotter space. 

From Eq. ([6|), one can define a local transfer operator 
r whose matrix elements are given by 



1,1+1 



>y, y +i - K 1 " < +1 |e- £h Hl - (7) 

An important step in the calculation below is to de- 
fine this local transfer matrix using fermion operators. 
Through a tedious calculation, we find that this transfer 
matrix can be expressed as an exponent of a quadratic 
function of fermion operators. 



1,1+1 
r y,y+i 



u k , y exp 



p k , y d\di +1 + qk,ydj +1 di 



-rk 



y (4 d l 



4 +1 d l+ i 



(8) 



where cf s are fermion operators defined in the Trotter 
space and coefficients (p k , y ,qk.y,i~k,y) are determined by 
the following equations 



sinhs^y 

Sk,y 

sinh s k y 

s k,y 
sinhsfc.j, 

'^k,y 
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Ik,- 



>'k. 



lk, y exp(q fc ^) 
eisinl^^ 

lk,y exp (-Qfc.y) 

et sinh7 fc ^ 

a k,y 
St ' 



Since coefficients (pk,y,q_k,y,rk, y ) do not depend on I, 
the above transfer operators are translationally invariant 
in every two unit cells along the Trotter direction. Thus 
one can block diagonalize these transfer matrices by tak- 
ing the Fourier transformation of fermion operators in 
the Trotter space. In the thermodynamic limit, Z k can 
then be simplified as 



Ny/2 



Z k = lim JjTr \{ (t k , 



2y-l, 



jtk,2y) ■ 



(10) 



where w = (2m+l)ir / M (m = 1, • • • , M) is the imaginary 
frequency. In obtaining the above expression, Fursten- 
berg's theorem [13] was used. t k ^ y -i,^ and t k .2 y are 2x2 
matrices defined by 



tk,2y-l,uj 
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with 



l k,y 



b k,y 



jk, y cosh ± a k ^ y sinhTfc^ 
-et sinh7fc, a 

lk,y 



= exp(±a kt y) 
-et sinh7 fe!a 



Thus the partition function Z k can be expressed as a 
product of N y 2x2 matrices. This is a great simplifica- 
tion to the problem, since the computer time needed for 
evaluating this product of 2 x 2 matrices scales just lin- 
early with N y . Furthermore, there is no need to store all 
these transfer matrices in advance. The computer mem- 
ory needed in the calculation is very small. Thus a truly 
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FIG. 1: Magnetization of the Hofstadter model at half-filling. 
The inset shows more clearly the lattice correction to the 
dHvA oscillation in the high field regime. 

big system with N y ~ 10 8 can be handled without any 
technical obstruct. 

In the Landau gauge, the lattice rotational symmetry 
is broken. The finite size effect along the x direction is 
small. In the temperature range we have considered, we 
find that N x = 50 is large enough. However, along the 
y-axis, the finite size effect is strong. We have evaluated 
the magnetization at T — 0.02 by varying N y from 500 
to 160000. We found that the results converge only after 
N y is above 50000. It indicates that indeed large lattice 
systems are needed in order to explore thermodynamic 
properties of the Hofstadter model. For higher temper- 
ature, the convergence can be reached with smaller N y . 
For the results shown in Figs. Q] and [21 we take N x = 50 
and N y = 80000 to ensure convergence. For simplicity, 
here we consider the half-filling case only. At half fill- 
ing, the chemical potential is pinned to /i = because 
of particle-hole symmetry. In the discussion below, the 
hopping constant t is set to 1 and e = 0.02. 

Fig. [1] shows the quantum oscillation of magnetization 
at three different temperatures. In the low field limit, the 
conventional dHvA oscillation is observed. The period of 
the oscillation, A(l/0), is about 2, consistent with that 
obtained from the textbook formula j 1 51] 

4tt 2 

A(l/<£) = — , (11) 

where Sf is the Fermi volume. At half filling, Sf = 
47r 2 /2 = 2-7T 2 . However, with increasing (j>, some sub- 
tle structures appear above the conventional dHvA curve 
expected for two dimensional electron gas in a magnetic 
field (see, the inset of Fig. [I}. These subtle structures 



become more and more pronounced with decreasing tem- 
perature. They result from the lattice correction to the 
energy spectra. 

Thermal fluctuation affects strongly on the line shape 
of magnetization. At high temperature, say T — 0.1, the 
fine fractal structure of Hofstadter butterfly with energy 
scale less than ksT is smeared out by thermal fluctua- 
tion. Only the conventional dHvA oscillation survives, 
except in the high field limit. However, at low tempera- 
ture, say T = 0.02, the fine structures of Hofstadter but- 
terfly with energy scales comparable to ksT will begin to 
influence the magnetic response of the system. It yields 
the sharp peaks or deeps observed in the magnetization 
curve in high fields. By further reducing temperature, we 
found that more and more peaks and deeps, even in the 
low field range, will emerge from the dHvA background. 

Around each sharp peak or deep, there is a change be- 
tween diamagnetism and paramagnetism with increasing 
temperature. For example, around (f> ~ 0.3, the magne- 
tization decreases with increasing (f> at T = 0.02 and the 
system is diamagnetic; whereas at T — 0.1, the magneti- 
zation increases with <f> and the system is paramagnetic. 
This change from para- to dia-magnetism is apparently 
due to the change of energy resolution since the energy 
spectrum is unchanged. It is a manifestation of the frac- 
tal structure of Hofstadter butterfly. 

Fig. [2] shows the field dependence of magnetic suscep- 
tibility x f° r between 0.27 and 0.4 at T = 0.1 and 
T = 0.02. At high temperature, T = 0.1, \ is para- 
magnetic. However, at low temperature, T — 0.02, \ 
oscillates strongly with <fi. It shows a series of local max- 
ima and minima, at which \ is positive (paramagnetic) 
and negative (diamagnetic) , respectively. These extrema 
appear when the magnetic flux takes some rational val- 





5/17 

A 


4/13 

h 


6 

1 


19 

6/17 

I 


4/11 

A 


5/13 

A 




3/10 


i 


16 


7/20 

7/22 

5 


i v 

/14 v 
14 3 

T= 

— T= 


1 

8 

0.1 
0.02 



FIG. 2: Magnetic susceptibility for the Hofstadter model at 
half filling. The values of (j> corresponding to local maxima 
and minima in the T = 0.02 curve are marked. 
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ues (j> = p/q (see the rational numbers given in Fig. [2j. 
The maxima and minima appear when q is odd or even, 
respectively. 

The appearance of these extremes results clearly from 
the interplay between periodic potential and magnetic 
field. It is strongly correlated with the density of states 
of the system at the Fermi level. At half filling, the den- 
sity of states vanishes linearly at the Fermi level (namely 
at a Dirac point) when q is evenQ. However, there is a 
van Hove singularity at the Fermi level and the density 
of states is divergent when q is odd[§|. Thus the mag- 
netic response is paramagnetic if (f> is close to a van- Hove 
singularity, and diamagnetic if <f> is close to a Dirac point. 

However, this connection between the extrema and 
density of states seems fragile if considering that there 
are infinite rational numbers p/q with even and odd de- 
nominators in an arbitrary small but finite interval of <f>. 
In other words, near any rational number, say <f> = 4/13, 
there are infinite other cf> at which the density of states at 
the Fermi level can be either zero or divergent. So how 
can we attribute the orbital paramagnetism at 4/13 to 
the van-Hove singularity in the density of states? 

This problem can be resolved by considering the hier- 
archical structure of Hofstadter butterfly and the temper- 
ature smearing of the band structure. At the first rank 
of hierarchy, the Hofstadter butterfly is divided into sev- 
eral subcellsJBj]. These subcells can be further divided 
recursively into many sub-subcells. This hierarchical re- 
cursion defines a parallel iterative transformation. After 
this iterative transformation, any rational <f> can be finally 
reduced to a simple rational number, which is equal to 
either 1/q' or 1 — l/<?', where q' is an integer. 

For example, <fi = 4/13 can be reduced to 4/5 after 
only one iteration. This means in the first order subcell 
centered at the Fermi level, there are five subbands and 
the middle one has a divergent density of states crossing 
the Fermi level. On the other hand, <j> — 401/1300, which 
is a value very close to 4/13, can be reduced to 3/4 after 
19 iterations. This means that, in the 19th order subcell, 
there are four subbands and the middle two meet at the 
Fermi level. In this case, the Fermi level is a Dirac point 
and the corresponding density of states vanishes. How- 
ever, the characteristic energy scale of the 19th order 
subcell is very small compared with the thermal energy 
ksT atT — 0.02. Therefore, the contribution by the sin- 
gularity at <j> = 401/1300 to x is completely smeared out 
by thermal fluctuation and only the peak at <f> = 4/13 
can be seen at T = 0.02. 

This can be seen more clearly by integrating out the 
density of states in an interval of fceT around the Fermi 
level. We find that the integral at <j> — 401/1300 is hardly 
different from that at <f> — 4/13. Therefore, around 
(j> ~ 4/13 the density of state at the Fermi level is de- 
termined by the van-Hove singularity in the first subcell 
at 4/13. Similar argument can be applied to 4> with even 
q. The difference is that in that case the density of states 



is dominated by the Dirac points. 

The above argument implies that the higher order sub- 
cells of Hofstadter butterfly can be probed by increasing 
the energy resolution. Thus more van-Hove singularities 
and Dirac points can be discerned by lowering temper- 
ature. However, in the limit of zero temperature, the 
susceptibility is no longer a well defined quantity, since it 
can oscillate between paramagnetism and diamagnetism 
in an infinitesimally small interval of 4>. This is consis- 
tent with the fact that that ground state energy is not 
differentiable with respect to <j>. 

In conclusion, we have introduced a quantum transfer 
matrix method to study thermodynamic properties of the 
Hofstadter model on square lattices. In the framework 
of Landau gauge, the partition function of the system 
can be determined simply by evaluating a product of a 
series of 2 x 2 local transfer matrices. This method al- 
lows thermodynamic quantities to be accurately and ef- 
ficiently evaluated without suffering from the finite size 
effect. It reveals many interesting features of quantum os- 
cillation in the magnetization. At high temperature, the 
lattice correction is small. However, at high field and low 
temperature, the lattice correction becomes significant 
and the magnetization shows many additional paramag- 
netic peaks or diamagnetic deeps above the background 
of dHvA oscillation. These peaks and deeps result from 
the fractal band structure of Hofstadter butterfly and 
are strongly correlated with the singular behaviors of the 
density of states at the Fermi level. 
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